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Abstract: We study Freeze Out process in high energy heavy ion reaction. The description of the 
process is based on the Boltzmann Transport Equation (BTE). We point out the basic limitations 
of the BTE approach and introduce Modified BTE. The Freeze Out dynamics is presented in the 
4- dimensional space-time in a layer of finite thickness, and we employ Modified BTE for the realistic 
Freeze Out description. 



I. INTRODUCTION 

The Freeze Out (FO) is an important phase of dy- 
namical reactions. It is of primary importance in case of 
rapid, dynamical processes where the originally strongly- 
interacting and locally equilibrated matter undergoes 
a rapid explosive process, in which matter properties 
change considerably, the interaction vanishes in a rela- 
tively small space-time layer, and local equilibrium dis- 
appears. The connection of the kinetic description of this 
process and the Boltzmann Transport Equation (BTE) 
raised considerable attention recently pjQ]- 

The problem is to calculate the phase-space (PS) dis- 
tribution of the post FO particles. Earlier such kinetic 
FO calculations were performed in one-dimensional mod- 
els tiylOjLij where the dynamics was governed by two 
constants: a re-thermalization parameter and a FO pa- 
rameter. This latter one is governed by the phase-space 
FO probability, which was constructed recently in a fully 
covariant form pj. 

The FO is a kinetic process and one would think it 
can be handled perfectly by using the Boltzmann Trans- 
port Equation, which may describe equilibrium and non- 
equilibrium processes equally well in a 4-dimcnsional 
space-time volume element, which is usually a FO layer. 
This work and ref. [2( follows this approach. This finite 
layer is frequently idealized as a 3-dimensional FO hy- 
persurface. In ref. 0,^3 author analyzes the features of 
this idealized discontinuity |27j . 

The FO can also be simultaneous with a phase tran- 
sition, especially when the phase transition reduces the 
number of degrees of freedom and contributes to the FO 
process this way. As an example let us describe a grad- 
ual hadronization and FO of the Quark-Gluon Plasma in 
a layer, where quasi-hadrons or hadrons are formed, the 
new particles gradually cease to interact, their PS distri- 
bution changes and the matter gradually freezes out. 

Free hadrons, which are formed, do not interact with 
anything and propagate directly to the detector. Al- 



though, the formation of these fragments can be most 
suitably described in a coalescence or recombination 
model, most finally observed baryon abundances follow 
the statistical model predictions. The reason is simple: 
the formation cross sections are governed by the same 
statistical factors as the thermal equilibrium, because 
the radial part of the formation probability for s-wave 
hadrons is about the same. Exceptions are the excited 
states, e.g. the p-wave hadrons like A(1520), which have 
a smaller radial form factor and consequently they are 
suppressed in comparison to the statistical equilibrium 
abundance, which is sensitive to their weight only. 

In this work we do not discuss issues related to simul- 
taneous hadronization and FO, for simplicity we con- 
sider one type of particles only and study their ki- 
netic evolution. If hadronization happens simultaneously 
with kinetic FO the kinetic description presented here 
can and should be extended. For example, some fea- 
tures of the fast hadronization and FO of supercooled 
Quark-Gluon Plasma, which might be created in ultra- 
relativistic heavy ion collisions^ for the first time such a 
scenario was proposed in refs. |3,|!j) are discussed in refs. 
[Tfl Hll | . The simultaneous FO and hadronization can be 
described in an idealized way by using the 3-dimensional 
FO hypersurface approach as suggested in ref. pj . Then 
this simplified approach makes it possible to solve some of 
the basic problems. The simultaneous hadronization and 
FO can also be handled by assuming idealized hadroniza- 
tion at the hypersurface, which is the inside boundary of 
the FO layer of finite thickness, L. In this way the present 
work is also relevant for the simultaneous hadronization 
and FO problem. 

In the present work we analyze the situation, discuss 
the applicability of BTE, and point out the physical 
causes which, limit the applicability of the BTE for de- 
scribing FO. And the aim is to show how can we overcome 
this obstacle. For this purpose we will modify the BTE 
and then will show how one can derive out of it a simple 
one dimensional kinetic model, similar to the one used 
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by some of the authors in earlier calculations. 



II. PARTICLES EMERGING FROM FREEZE 
OUT HYPERSURFACE 

Not only in heavy ion reactions, but in many dynami- 
cal processes particle creation (or condensation) happens 
mostly in a directed way: the phenomenon propagates 
into some direction, i.e. it happens in some layer or 
front (like detonations, deflagrations, shocks, condensa- 
tion waves or FO across a layer with space- like normal). 
The reason is that neighboring regions in the front may 
interact to minimize the energy of the front by evening 
it out, providing energy to neighboring regions to ex- 
ceed the threshold conditions. Even in those relativistic 
processes that are time-like (have time-like normal), and 
so the neighboring points of a front cannot be in causal 
connection, the dynamical processes may and frequently 
have a direction. See the example in ref. 13] . This 
can be a simple consequence of the initial and boundary 
conditions. 

These fronts have a characteristic direction (or normal, 
da' 1 ). Let us look at an example when particles in a 
domain of the space-time (ST) are characterized by a 
phase-space distribution, f(x,p). Then the space-time 
current density of these particles, N ,1 (x) can be described 
as 

V" / fPpl>f(x,p) . (1) 



P 

The net number of particles crossing an arbitrary hyper- 

d 3 p 



surface element da^ is 



dS = N^da^ = 



P^ da^ f(x,p) 



(2) 



If we want to describe the FO, particles are allowed to 
cross the FO hypersurface "outwards" only, i.e., only in 
the direction of da^ . Thus, 



~>FO 



N Fo da t* 



(3) 



d 3 p 
IF 



p** da^ fFo(x,p)Q(p' J ' da^) 



where either the phase space distribution, fFo(x,p), 
should have only particles with momenta pointing out- 
wards (post FO distribution), and/or this is secured by 
the step function Q(p fJl da^). Eq. (@J yields the mod- 
ified Cooper- Frye FO formula, where Jfo(x,p) should 
be determined in such a way that all conservation laws 
across the FO hypersurface are satisfied and overall en- 
tropy does not decrease! 0, Q, 



III. NON-ISOTROPIC PARTICLE SOURCES 



and it is more realistic to assume a continuous, 4- volume 
FO in a layer (or domain) of the space-time. At the in- 
side boundary of this layer no particles are frozen out yet, 
while at the outside boundary hypersurface all particles 
are frozen out and no interacting particles remain (see 
Figure For the sake of simplicity let us also assume 
that the total particle number is conserved, even if simul- 
taneous freeze out, hadronization and particle formation 
are frequently discussed. 

Thus, while the total number of particles remains con- 
stant, in this domain, the number of interacting particles 
decreases and the number of frozen out or free particles 
increases: 



N»(x) 
dpN»(x) 



d,K(x) 







(4) 
(5) 
(6) 



Then the space-time (ST) volume element, d 4 x, in the 
layer of interest can be converted into d 4 x — > ds^ da^, 
where ds 11 is the length element in the direction of the 
4- vector da 11 , which can be space- like or time-like, i.e.: 
time-like, da 11 da^ = +1, or space-like, da' 1 da^ = — 1. 
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FIG. 1: Space-time picture of the FO process. At early 
times, centrally in the collision region we have intensively in- 
teracting matter, which is equilibrated and thermalized, this 
is the fluid-dynamical domain bordered by the S2 hypersur- 
face, which has a normal 4- vector dcr M . The location of this 
surface is given by the fact that the interacting fluid is cooling 
and expanding, and reaches a point when interactions are not 
frequent enough to maintain full thermal and hydrodynami- 
cal equilibrium locally. Some particles will not interact any 
more beyond this hypersurface. Later on in the expansion 
and cooling we reach an other hypersurface, Si. By reaching 
this surface on their way all particles become non-interacting, 
or free. Thus, when reaching this surface the FO process is 
completed. The momentum distribution of particles does not 
change any more. This is the (post) FO distribution 



The FO-fronts or FO-layers are not necessarily narrow, 
but they have a characteristic direction (or normal, da 11 ), 



Let us also assume that the ST domain, where free 
particle formation happens, is a layer, which is relatively 
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narrow compared to the bulk of matter (see Figure^. 
Assume also that the boundaries of this layer are parallel 
or approximately parallel, and the thickness of the layer 
does not vary much. Under these conditions one can 
describe the change of free particle number in the layer 
via the divergence of the 4-current of the particles by the 
expression: 



AJY; 



d 4 x dpNfix) 



(7) 



= / ds 11 da,, 



d 3 P ~o d t*fi( X 'P) ■ 



According to the physical assumptions discussed above, 
the 4-divergence is maximal in the direction of da v , and 
negligible in the other 3 orthogonal directions. 

The emission or freeze out probability may depend 
on physical processes, cross sections, transition rates, 
and the actual PS distributions. Furthermore, fi can 
be space-time dependent, and must be determined self- 
consistently durin g th e detonation, deflagration or FO 
process HEHHlll. 

We will return to realistic FO probabilities later in sec- 
tion |W1 



IV. FREEZE OUT AND THE BOLTZMANN 
TRANSPORT EQUATION 

One can derive the Boltzmann Transport Equation 
from the conservation of charges in a ST domain [l7| . 
A 4 x, assuming the standard conditions: (i) only bi- 
nary collisions are considered, (ii) we assume "molecular 
chaos" , i.e. that the number of binary collisions at posi- 
tion x is proportional to f(x,pi) x f(x,p2), and that (hi) 
f(x,p) is a smoothly varying function compared to the 
mean free path (m.f.p.). The conservation laws lead then 
to the requirement that the integral of the 4-divergence 
of conserved charges should vanish 



f f d 4 x^-p»d fl f(x ) p) = 0. 
JA*x J ah P 



(8) 



As the choice of the ST 4- volume element is arbitrary we 
obtain the differential form of the conservation law, which 
describes the evolution of the PS distribution, f(x,p), of 
a particle with momentum p. However, if we take into 
account that particles can scatter into this PS volume 
element around p, or can scatter out from this volume 
element, we have to add Gain- and Loss- collision terms 
to the conservation equation (see e.g. sect 3.2 of ref. 

E2): 



i 2 23 4 f( Pl )f(p 2 )W^ 2 
2V 3 J(p)f( P2 )W^ . (9) 



Here we assume elementary collisions where in the initial 
state two particles collide with momenta p\ and p 2 into 



a hnal state of two particles with momenta p 3 and p±. In 
case of the Gain term the particle described by the BTE, 
with momentum p (without an index), is one from the 
two final state particles, while in case of the Loss term 
this particle is one of the initial state particles. This is 
indicated by the indexes of the invariant transition rate 
[l7j . We integrate over the momenta of the other three 
particles participating in this binary collision. We use 
the notation 
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2? a = 



d 3 pi d 3 p 2 d 3 p 3 



P\ 



P% 



We can shorten the notation further by suppressing the 
arguments of the PS distribution functions, and the in- 
dexes of the momenta in the argument will be carried 
by the distribution function / and the collision term 



P1P2 



TP^dJ 



Now, aiming for the description of FO process let us split 
up the distribution function, /, into / = /' + /* , where 
ft is the phase-space distribution function of the "free" 
or frozen out particles, which do not collide any more, 
while f l is the interacting component 0,13. Thus, the 
FO process is represented here by gradually populating 
and building up the "free" component, while draining 
particles from the interacting component. As the parti- 
cles belonging to the free component may not collide any 
more, they do not appear in the initial state components 
of collision integrals! 



P^if + ft) 



The gain term, f\ f\ W± 2 populates both the interact- 
ing, f l , and free, /*, components, so we will introduce 
a FO probability, which 'feeds' the free component. The 
probability is phase-space dependent. In principle it may 
depend on the positions and momenta of both incoming 
particles, and it can weight the outgoing phase space for 
one (or both) outgoing particles. In the most simple case 
we have to assume that it depends at least on the mo- 
mentum of the outgoing particle, which belongs to the 
component /* : V F (x,p) = V$. 
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W + ft) = - I uVJlfl V f W£ + (1 - V f )W, 



rpi 



rpi 



2 2W72%2 4 - 



(12) 



Now, we can separate the two components into two equa- 
tions. The sum of these two equations returns the com- 
plete BTE above: 



(13) 
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(14) 

The free component does not have a Loss term, because 
particles in the free component cannot collide, and so, the 
free component cannot loose particles due to collisions. 
Rewriting the second equation yields: 

P^f = -\j i2V A fir 2 V f W^+ l - J uDtfifiWg 

- \ J 2D M ffiW% (15) 

The first term is a drain term, describing the " escape" or 
"freeze out" of particles from the interacting component. 
It is the inverse of the gain term (or source term) for 
the free component, j' . The last two terms are influenc- 
ing the interacting term by redistributing particles in the 
momentum space. These latter two terms do not include 
the FO probability factors! Thus, these two terms drive 
the interacting component towards re-thermalization. As 
a usual approximation these two terms can be approxi- 
mated by the relaxation time approximation as in refs. 
HHEl. Thus ' the BTE describing FO in this situation 
reads as: 

P^J f = \J vfflififtPjWg (16) 
P^f =~l i*D A fiftp f Wg + p° (17) 

I J T re l 

The first equation, eq. ifH^I. describes the gain of the 
free component, i.e. that part of the earlier gain term, 
which will not collide any more. The first term in the 
second equation has the same value with opposite sign. 
This describes the part of which is leaving the in- 
teracting component and does not take part in the re- 
thermalization. 

In fact the above described collision integrals can be 
highly simplified, by exploiting the symmetries and con- 
servation laws in the invariant transition rate, W, so that 
only one phase-space integral remains to be executed (see 
section 3.3 and eq. (3.27) in ref. 01)- 

V. MODIFIED BOLTZMANN TRANSPORT 
EQUATION 

Now, the question arises: can the BTE handle real- 
istically the FO process? We have seen that the struc- 
ture of the kinetic equations, used earlier to describe FO 
0,0,0,0, and the separation of the "escape" and "re- 
thermalization" terms come out in a simple, straightfor- 
ward way from the BTE. 

However, the usual structure of the collision terms in 
the BTE are not adequate for describing rapid FO, in 
a layer which is comparable to the m.f.p. If we assume 
the existence of such a layer this immediately contradicts 



assumption (iii): the change is not negligible in the di- 
rection of do v . The assumption (ii) of "molecular chaos" 
is also violated in a FO process because number of colli- 
sions is not proportional with f(x,pi) x f(x,p2), but it 
is delocalized in the normal direction with f(x + \,pi) x 
fix — X,P2)- (The fact that the FO is a delocalized ki- 
netic process, was already used in ref. Q when integrals 
along the path of propagating particles were introduced, 
but the consequences regarding the details of the collision 
terms and the validity of the molecular chaos assumption 
were not discussed.) 

Based on the above considerations, one might conclude 
that the changes of the distribution function are medi- 
ated by the transfer of particles, and consequently only 
slowly propagating changes are possible. I.e., the front 
propagates slowly, and its normal, da 1 *, is always space- 
like. This was a common misconception, until recently, 
where all " superluminous" shock, detonation, deflagra- 
tion fronts or discontinuities were considered unphysical 
based on early studies |l8| . However, it was shown re- 
cently, that discontinuous changes may happen simulta- 
neously in spatially neighboring points, i.e. the normal of 
the discontinuity- hypersurface can be time-like |l3l |19| . 
This applies to the FO process also. Thus, the direction 
of characteristic or dominant change, da 11 , may be both 
space-like and time-like in the FO process. 

From all the processes mentioned above (i.e. shocks, 
detonations, deflagrations etc.) the FO is the most spe- 
cial one. Because the number of interacting particles 
is constantly decreasing as the FO proceeds and corre- 
spondingly the m.f.p. is increasing and, in fact, it reaches 
infinity when the complete FO is finished. This simply 
means that we strictly speaking can not make FO in fi- 
nite layer of any thickness smooth enough to be modeled 
with BTE. It is also obvious that if FO has some charac- 
teristic length scale (thickness of the layer or even some 
characteristic parameter for infinitely long FO 4]), it is 
not proportional with the m.f.p., because m.f.p. increases 
as the density of interacting component becomes smaller, 
while FO becomes faster in this limit, so its characteristic 
scale should decrease. 

Since, there is a strong gradient in the FO direction: 
the free component rapidly increases, while the interact- 
ing component decreases along the FO direction, we can 
conclude that the collision terms in their usual form are 
not adequate to describe the FO process, particularly not 
the "escape" probability or " escape" term. The appropri- 
ate equations to describe this system can be a Modified 
Boltzmann Transport Equation (MBTE) [H : 

p'd^fip) = \J uVjix^T) x f&plj X?p 2 

- \ J iD^JM) X W~P2) (18) 



where f(x,Pi) is an average over all possible origins of 
the particle in the backward lightcone of the ST point 
x = (t, x): 
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/* dii/d 3 x 1( 5 3 (f - xt - v(t - t 1 ))f(x 1 ,p)e- f n *»/^*W* a )«« 8 (*»-ft-^* a -ti)) 
// dt 1 /d 3 x 1 5 3 (f-f 1 -^(i-f 1 ))e-^ dt2/d3x2CT " (a2) ^ 3(S2 ^ 1 ^ (i2 - tl)) 



(19) 



where <5 3 (x — xi — u(t — ti)) fixes the ST trajectory, along 
which the particles with given momentum can reach the 
ST point x, time to is given by the initial or boundary 
conditions, v = p/p (v = \v\), and the exponential fac- 
tor accounts for the probability not to have any other 
collision from the origin x± till x. In the arguments of 
exponents n(x) is the particle density in the calculational 
frame, n(x) = N°(x), and a is the total scattering cross 
section. After performing integrations over d 3 x with a 
help of (5-functions we can write the MBTE equation in 
the form: 
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V t / 2 f(t 1 ,p 1 )G(t uPl )f(t 2 , P2 )G(t 2 ,p 2 )WPP p2 



2 / 2VlTf(tuP)G(h,p)f(h,P2)G(h,p 2 )W^ 



where 



12^4 2 = o / d *l / *a / 12^4 , 



/(*1,P) = f(ti t x-v(t - ti),p) , 



(20) 



(21) 



(22) 



— J t dt2<Tn(t2,x—v(t — t2))v 



C{x,p) = ( dt 1 e~ I *i dtatm{t " g ~^ t - t%))v . 



(24) 



to 



Interestingly, Molecular Dynamics models do not use 
the local molecular chaos assumption, and follow the tra- 
jectories of the colliding particles instead. Thus such 
models do actually solve the MBTE, and not BTE, al- 
though this was not realized before. 

The obvious limit in which MBTE is reduced to BTE is 
a completely homogeneous ST distribution function (i.e. 
no external forces, no boundaries). Another possibility 
is the hydrodynamic limit, A = 1/an — > 0, when the 
exponential factors H23I24|> will be reduced to ~ S(t — 
£1,2), reproducing the BTE after t\,t 2 integrations. 

The symmetries and the assumption of local molec- 
ular chaos lead to the consequence that local conser- 
vation laws can be derived from the original BTE, i.e. 
d^T^ = and d^N 1 * = 0, where T and N are given 
as integrals over the single particle PS distribution and 



the momentum. Although now we have delocalized the 
equations, the local conservation laws can still be de- 
rived in the same way, as it was shown in Ref . |20| . The 
very essential property of the BTE is the Boltzmann H- 
theorem. Here the situation is more complicated and the 
behaviour of the entropy current in MBTE is a subject 
of future studies. Nevertheless, for adiabatic expansion 
(Sjt = 0) a sufficient condition is the same as for BTE, 
namely f(x,px)f(x,p 2 ) = f(x,p 3 )f(x,p 4 ) H3- 

The obtained MBTE is considerably more compli- 
cated, than the original BTE. In order to proceed let 
us make a further simplification assuming that all the 
particles arrive into the collision point x from one m.f.p. 
distance |2l| (instead of allowing them to arrive from 
any distance with the corresponding probabilities, as it 
is done in eq. (119(1 ). This then leads to the following 
simplified equation, which nevertheless is still adopted to 
the strongly non-homogeneous systems much better then 
the original BTE: 



P^fip) = \( 12 VJ(x 1 ,p 1 )f(x 2 ,p 2 )W« 



PP4. 



iV 3i f(i,p)f(x 2 ,p 2 )W^, (25) 



where Xk is the origin of colliding particles, i.e. the ST 
point where the colliding particles were colliding last, 
Xk = x - u k T(x,Vk/vk), Ufc = (7fe,7feWfc), 7 = l/Vl - v 2 
and Vk = Pk/Pk- Here r is the collision time, such that 
\v\t{x,v /v) — X(x,v/v). Note that the m.f.p. depends 
not only on the position, but also on the direction of the 
particle motion. This is an essential modification if the 
PS distribution has a large gradient in the space-time. 
This gradient defines a ST 4-vector characterizing the 
direction of the process, dcr^. In ref. the direction 
du^ is also introduced, however, it is not discussed why 
and it is not connected to the derealization of the BTE. 

For the FO modeling, repeating for the eq. I|25() the 
same step as for BTE above, we then obtain: 

p^ff(x,p) = i J 12 V 4 V f W^ f i (x 1 ,p 1 )f(x 2 ,p 2 ) , 

(26) 

P^dJ\x, P ) = -\J nDPfWg f (xi^O/* (ai 2 ,p 2 ) 

->eq J 



+ p" 



(27) 



A simple general solution of the MBTE I|2()I27|) can- 
not be given but it serves as a basis for simplified, phe- 
nomenological kinetic models describing the FO process. 
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VI. APPROXIMATE KINETIC FREEZE OUT 
MODELS 

In this section our goal is to present a schematic deriva- 
tion of a simple kinetic FO model used by some of the 
authors earlier. This represents only one particular pos- 
sibility and the general MBTE equation can be solved 
or approximated in other ways also. The approximation 
we present is one of the simplest possibilities, but not 
necessarily the most realistic one. 

If the ST distribution is non-uniform and the direction 
of steepest gradient can be clearly identified, one may 
replace one (or more) of the integrals over d 3 pi (or d 3 p2) 
by space-time integrals over the origins of the incoming 
particle(s), d 4 x', requiring that the particle reaches the 
ST point, x, when needed. This requirement determines 
p M for a given x /fl . It is reasonable to assume that after 
converting some of the integrals to ST integrals and per- 
forming them, we get an effective FO term reflecting the 
properties of the local PS distribution, transition rate, 
the ST configuration (e.g. gradient of density change, 
and its direction) and characteristics of the FO layer. 

Let us return to the basic integral form of the kinetic 
theory, eqs. (|26I27[) . and discuss the FO probability We 
will study equation i|2rj|) without performing the integrals 
in a formal way, rather illustrating the procedure giving 
a better insight into the problem. 

When we are in the FO layer, close to the bound- 
ary of complete FO we have to calculate here the col- 
lision rate. According to the MBTE this depends on 
the PS distribution of the incoming particles at their ori- 
gins, f l (xi,pi) f l (x2,Pi)- Assume that the FO direction 
points in the direction of da 11 , as it is shown in Figure 
13 On the right hand side of the collision point the den- 
sity of interacting particles is low or zero, while on the 
left hand side it is larger, closer to the pre FO value (see 
Figure 01 . 



It is more probable, that particles arrive to the collision 
point x from the left side, because of the higher density 
of the interacting particles on the left. Consequently, 
most outgoing particles leave to the right. Thus, the col- 
lision rate at x depends on the conditions what we have 
around i\ and X2, i.e. deeper inside the interactive mat- 
ter. Consequently, the collision rate is still higher than 
the conditions at x could secure! Then V f determines 
what fraction of the outgoing particles will freeze out 
from those which collided. The collision rate does not 
go to zero even if we are at the outside boundary of the 
FO layer, because particles still can arrive from the left 
where we still have interacting particles. As there are no 
interacting particles on the right hand side, all of these 
particles should freeze out, i.e. Vf — > 1, when x — > L 
(see Figure |2Jl- 
Let us execute two of the phase space integrals for one 




FIG. 2: The plot of one of the last collisions at x, shown 
in the spatial cross section of the FO layer. Particles arrive 
from positions xi and xi to point x with momenta pi and 
P2- Within the FO layer of thickness, L, the density of in- 
teracting particles gradually decreases (indicated by shading) 
and disappears at the outside boundary, Si (thin line) of the 
layer. R.h.s. from this boundary there are no interacting 
particles. Particles can reach x from a region closer than the 
mean free path (m.f.p. indicated by the dashed line), but only 
from places where the interacting particle density is still not 
zero, i.e. mostly from the left. The inside boundary of the FO 
layer, S2 (thick line) indicates the points where the FO starts. 
Left of this line there is only interacting matter and the FO 
probability is assumed to be zero for collisions happening in 
the interacting region. 

incoming and one outgoing particle, J ^tr- 

pfdpf* = \J nVJin V f W$ = (28) 

= \QM J ^rfPuPi) PfW*, (29) 

where Q2 = J f l {xi,P2) and V4 are invariant scalars. 
Eq. J2HI resembles eq. (3.27) in ref. 17], but one of the 
incoming particle distributions, f l {x2,P2) is integrated 
out, and leads to an integral quantity, Q%- This can be 
approximated by the invariant scalar density at X2, i.e. 

Q 2 « n 2 (x 2 ) . 

Here V4 is not known directly, but can in principle 
be calculated based on the distributions, f % (x\,pi) and 
P(x2,P2), and the transition rate, W. The resulting 
transition rate will then be averaged over particles 2 and 
4 W pi - 

p»d fl ff(x,p) = lQ 2 V i J ^fixt^VfWg. (30) 

As we see this resulting equation is delocalized - in a 
rapid dynamical process the distributions at x and X\ 
are not the same, as it was discussed above. 
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Now, eq. (|30fl can be integrated either in the x\ - 
space, or in the p\ -space, as the two are connected by 
the fact that a particle should travel from x% to x with 
momentum p\. We should integrate over all x\ points 
from where one can reach x in a collision time, t(x, v/v). 
This brings in information about the local spatial gra- 
dient of the ST distribution function, as we discussed it 
above. The collision terms in the original BTE contain 
only local information, which is assumed to be isotropic 
(or slowly changing), so it is neglected. 

In addition the FO probability, Vf, may include inte- 
grated information about the FO process, e.g. the prob- 
ability not-to-collide with anything on the way out, rea- 
sonably should depend on the integral number of inter- 
acting particles on the way out. 

For the sake of simplicity let us assume small an- 
gle scatterings, and the propagation of a single particle 
~ tt>| 5(p — pi), then 

jrd»ft{x,p) = f{x x ,p) V f w\ . (31) 

The cumulative effect of all particles which can reach the 
ST point x in a collision time, leads to a change directed 
into the direction given by da^. The transition rate, 
can be estimated as (av re i) ~ p^/p , what yields to: 

p»d„ff(x,p) = f(x uP ) ^Q»V A V f da^/p° .} 

(32) 

As we mentioned, the spatial variation of the phase- 
space distribution cannot be neglected in rapid dynami- 
cal processes as the FO, and this brings in a direction of 
the dominant change, da' 1 . 

Let us now consider the FO situation, where we have 
a directed process in a layer. The dominant change hap- 
pens in the direction of the normal of the FO hypersur- 
face, da' 1 (where da^ da^ = ±1). We can decompose the 
4-vector, p v on the l.h.s. of the above equations into four 
orthogonal directions: 

p v = (p" do ^)da v + (p" daedal (33) 
+ (p tl da2„)da^ + (p^ da 3 „)da% , 

where the 4- vectors, da\, da%, and da 3 , are tangent 
to the hypersurface and orthogonal to the normal, da v . 
This leads to: 

p v d v f{x,p) = [(pPfoJdo" + W&nJdoZ (34) 
+ (jfda 2li )da v 2 +(jfda 3lx )dal]d v f{x,p) . 

Here we assumed that the change happens in the direc- 
tion of the normal and negligible along the hypersurface 
of the front, thus the last three terms can be neglected: 

p"d v f(x,p) w {p^da^) da u d u f(x,p) . 

Inserting the above equation into llol't yields a kinetic 
equation describing the directional derivative of the 



distribution function in the direction of the dominant 
change, da^ as 

da»d^{x,p) = f\x 1: p)P: sc , (35) 

where the escape probability depends on the ST co- 
ordinates, on the interacting part of the PS distribu- 
tion, on the transfer properties and the FO probabil- 
ity: P* sc (x,p, /% da, w, Vf). The x% in this case means 
x\ = x — da 1 * Ap, where A p is a four- vector of m.f.p. 

This x\ in the argument of the distribution function 
on the r.h.s. of the eq. (|35|l is extremely important. 
Certanly, we can repeat all the steps from eq. l|28l) to 
eq. (|35|l based on the BTE amd the result will be the 
same except for the derealization of / on the r.h.s. This 
x—da^X^ dependence of the distribution function reflects 
the FO property which was discussed at the beginning 
of this section and illustrated in Fig. [21 - the collision 
rate at some point x, and correspondingly the number of 
particles, which will freeze out after this collision, feels 
the properties of the matter deeper (by about one m.f.p.) 
inside the interacting matter. 

The derivation above did neglect several details and 
features, however, reflects the basic structure of ad hoc 
kinetic FO models 0. 0. Q In these models the in- 
finitely long FO was studied, and therefore the dereal- 
ization of eq. I|35|l was not so important. For the FO 
modeling in the finite layer (2^ this effect will cause a 
substantial difference making FO faster. 

VII. ESCAPE PROBABILITY 

The escape probability in eq. (|35|l can be estimated 
based on fundamental physical principles, like it is done 
in the above mentioned works. The approach can, nev- 
ertheless, be improved if we take into consideration the 
origin of the above derivation, especially the requirement 
of full covariance of the model and the requirement that 
the FO process may point in any space-time direction. 
The first significant advances, where these principles were 
applied, are presented in |7J ■ It incorporates the achieve- 
ments of recent years, by cutting negative contributions 
in the FO density [l5j and making the FO direction de- 
pendent 0. Here we just present briefly a direct estimate 
for the escape probability [2^, |2|| . 

The escape probability includes the FO probability, 
which separates from among the outgoing (gain) parti- 
cles, which fraction of them is still colliding and which 
not. The probability not-to-collide with anything on the 
way out, reasonably should depend on the number of par- 
ticles, which are in the way of a particle moving outwards 
in the direction p/p, across a FO layer of estimated thick- 
ness L (representing the fact that we have finite number 
of particles on the way out to collide with 4] ) . If we are 
in this FO layer and progressed from the beginning of the 
layer to a position x^ , there is still 

L — x^da^ 
cos 6 



8 



distance ahead of us, where is an angle between the 
normal vector and p/p. We assume then that the FO 
probability is inversely proportional to some power of 
this quantity [HHIl. Thus 

where the power a is influencing the FO profile across 
the front, and the cut factor is eliminating negative con- 
tributions to FO. In papers [J [J 00 the authors have 
used o = l, and modeled FO in an infinite layer, i.e. in 
L — > 00 limit. Furthermore, they were using a constant 
characteristic length A instead of A(xi): 

C^ew. (37) 

Comparing eqs. l|3lfl) and (|3T|) one can see that now we 
replace the constant characteristic length A, which was 
clearly oversimplifying the situation, with two factors. 
The first is the collision rate, which is proportional with 
jrt-y « (n(xi)a) , and this does not tend to zero even 
if we reach the outside boundary of the FO layer, as 
this parameter is characteristic to the interior region at 
x\. The other is the generalized FO probability, which 
depends on the direction of the outgoing particle and 
on the number of interacting particles left in the way to 
collide with, i.e. oc y^— , where we have fixed da^ = 
(0, 1,0,0). So, we have generalized eq. lj57|) by replacing 

A — > X'(x) = \(xi)^-—^- . (38) 

Now the new characteristic length X'(x) gradually de- 
creases as FO proceeds and the number of interacting 
particles becomes smaller and smaller, and goes to 
when the FO is finished, as it was discussed in section 

El 

The simple angular factor, cosO, maximizes the FO 
probability of those particles, which propagate in the di- 
rection closest to the normal of the layer, da^. The quan- 
tities, cos O = p x /\p\ for FO in x-direction and cos = 1 
for FO in i-direction, are not Lorentz invariant. There- 
fore, to make our description completely invariant we 
shall generalize it to p » — ~ cos O. 

So, we write the invariant escape probability, within 
the FO layer cove ring both the time-like and space-like 
parts of the layer |2g , as 

(39) 

If we take the four velocity equal to u M = (1,0,0,0), 
in the Rest Frame of the Front (RFF), i.e. where 
da^ = (1, 0, 0, 0), then the momentum dependent part of 
the Escape Probability, P(p), is unity. Otherwise, in the 
Rest Frame of the Gas (RFG), where — (1, 0, 0, 0), the 
escape probability P(p) is P(p) = p^da^ Q^p^da^/p . 



More detailed investigation about escape probability 
P(p) for different da^ can be found in 0. 

To calculate the parameters of the normal vector da^ 
for different cases listed above, we simply make use of 
the Lorentz transformation. The normal vector of the 
time-like part of the FO hypersurface may be defined as 
the local i'-axis, while the normal vector for the space- 
like part may be defined as the local x'-axis. As the 
da^ normal vector is normalized to unity its components 
may be interpreted in terms of 7^ and v a , as da^ = 
7o-(l, Vcr, 0, 0), where j a = , 1 for time-like normals 

and 7 CT = , 1 for space-like normals. 

The detailed results of the application of this covariant 
escape probability will be presented elsewhere |22( • 

In refs. 0, || the post FO distribution was evalu- 
ated for space-like gradual FO in a kinetic model. Ini- 
tially we had an equilibrated, interacting PS distribu- 
tion, ji nt {p,x), and an escape probability, similar to eq. 
(|39|l . but simplified one. It was dependent on the an- 
gle of the two vectors only. After some small fraction of 
particles were frozen out as the FO process progressed in 
the front, the interacting component were re-equilibrated 
with smaller particle number, smaller energy and mo- 
mentum to account for the quantities carried away by 
the frozen out particles. This was then repeated many 
times in small steps along the FO front and the frozen 
out particles were accumulated in the post FO PS dis- 
tribution, ffree- The resulting distribution was highly 
anisotropic and obviously non-equilibrated. The details 
of the post FO distribution depend on the details of the 
escape probability, and on the level of re-equilibration of 
the remaining, interacting component. 

Bugaev assumed earlier (see also ref. [23|), that 
the post FO distribution is a (sharply cut) " Cut- Jiittner" 
distribution, but the above mentioned model shows that 
this can only be obtained if re-equilibration is not taking 
place. The kinetic model provided an asymmetric but 
smooth PS distribution while the escape probability 
(|39|l yields a somewhat different, but also smooth PS 
distribution |26| . These can be well approximated by the 
"Canceling Jiittner" distribution |24j . 

In ref. 6] the same infinite ID model (as in refs. 0,0) 
was applied for the FO through the layer with time-like 
normal vector. The model in this case can be solved ana- 
lytically (since cosO = 1) and thus, the exact form of the 
post FO distribution for pions and protons was obtained. 
Although analytical expressions for these distribution are 
different from the thermal Jiittner distribution, the forms 
of the functions are very similar for intermediate and high 
momenta. Deviations at the low momenta seems to be 
due to infinite long FO (they become much smaller for 
the escape probability l|39|l , modeling FO in a finite layer) 
and ID character of the model prif . 
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VIII. CONCLUSIONS 

The FO process was discussed in the 4-dimensional 
space-time in a layer of finite thickness. Arising from 
the physical process this layer is directed, it has an in- 
side and outside boundary, which are not identical. The 
processes in the layer are not isotropic, they must be sen- 
sitive to the direction of the layer. It is shown that, as 
a consequence, the basic assumptions of the Boltzmann 
Transport Equation are not satisfied in this layer, and the 
equation should be modified. It is also shown that earlier, 
ad hoc kinetic models of the FO process, can be obtained 
from this approach in a fully covariant way, and freeze 



out in space-like and time-like directions can be handled 
on the same covariant footing. 

Acknowledgments: One of the authors, L.P. Cs., 
thanks the Alexander von Humboldt Foundation for ex- 
tended support in continuation of his earlier Research 
Award. The authors thank the hospitality of the Frank- 
furt Institute of Advanced Studies and the Institute for 
Theoretical Physics of the University of Frankfurt, and 
the Gesellschaft fur Schwerionenforschung, where parts 
of this work were done. Enlightening discussions, with 
Cs. Anderlik, K.A. Bugaev, F. Grassi, Y. Hama, T. Ko- 
dama and D. Strottman are gratefully acknowledged. 



[1] K.A. Bugaev, nucl-th/0401060 

[2] Yu.M. Sinyukov, S.V. Akkelin, Y. Hama, Phys. Rev. Lett. 

89 (2002) 052301; S.V. Akkelin, M.S. Borysova, Yu.M. 

Sinyukov, [nucl-th/0403079] 
[3] L.P. Csernai, Zs. Lazar and D. Molnar, Acta 

Phys. Hung. New Ser. Heavy Ion Phys. 5 (1997) 467; Cs. 

Anderlik et al., Phys. Rev. C59 (1999) 388. 
[4] V.K. Magas et al., Acta Phys. Hung. New Ser. Heavy Ion 

Phys. 9 (1999) 193. 
[5] V.K. Magas et al, Nucl. Phys. A661 (1999) 596; Cs. 

Anderlik et al., Phys. Lett. B 459 (1999) 33. 
[6] V.K. Magas, A. Anderlik, Cs. Anderlik and L.P. Csernai, 

Eur. Phys. J. C30 (2003) 255-261. 
[7] L.P. Csernai et al., |hep-ph/0401005| 
[8] T. Csorgo and L.P. Csernai, Phys. Lett. B333 (1994) 494. 
[9] L.P. Csernai and I.N. Mishustin, Phys. Rev. Lett. 74 

(1995) 5005. 

[10] A. Keranen, L.P. Csernai, V. Magas and J. Manninen, 

Phys. Rev. C 67 (2003) 034905. 
[11] L.P. Csernai et al., Acta Phys. Hung. New Ser. Heavy Ion 

Phys. 17 (2003) 271. 
[12] K.A. Bugaev, Phys. Rev. Lett. 90 (2003) 252301. 
[13] L.P. Csernai, Sov. JETP 65 (1987) 216; Zh. Eksp. Theor. 

Fiz. 92 (1987) 379. 
[14] F. Cooper and G. Frye, Phys. Rev. D10 (1974) 186. 
[15] K.A. Bugaev, Nucl. Phys. A606 (1996) 559. 
[16] Cs. Anderlik et al., Phys. Rev. C59 (1999) 3309. 
[17] L.P. Csernai: Introduction to Relativistic Heavy Ion Col- 
lisions, Willey (1994). 
[18] AH. Taub, Phys. Rev. 74 (1948) 328. 
[19] K.A. Bugaev, M.I. Gorenstein, W. Greiner, Acta 

Phys. Hung. New Ser. Heavy Ion Phys. 10 (1999) 333-343; 

K.A. Bugaev, M.I. Gorenstein, W. Greiner, J. Phys. G25 

(1999) 2147. 

[20] . V.K. Magas, L.P. Csernai, E. Molnar, A. Nyiri, and K. 



Tamosiunas, Nucl. Phys. A749 (2005) 202. 
[21] L.P. Csernai, V.K. Magas, E. Molnar, A. Nyiri and K. 

Tamosiunas, hep-ph/0406082 

[22] E. Molnar et al. (2005) in prep aration. 

[23] K.A. Bugaev, M.I. Gorenstein, nucl-th/9903072 

[24] K. Tamosiunas and L.P. Csernai, Eur. Phys. J. A20 

(2004) 269-275. 

[25] V.K. Magas, Talk at the 4th Collaboration Meeting on 
Atomic and Subatomic Reaction Modeling + BCPL User 
Meeting, Trento, Italy, February 26-29, 2004. 

[26] V.K. Magas, Talk at the International Workshop "Cre- 
ation and Flow of Baryons in Hadronic and Nuclear Col- 
lisions" , Trento, Italy, May 3-7, 2004. 

[27] In refs. 0, ll2T | the author discusses two physically dif- 
ferent situations, which should be clearly separated, oth- 
erwise one can run into some confusions. The first one 
is the "transition" from hydrodynamical description of 
the system created in heavy ion collisions to cascade de- 
scription. This is not a physical phase transition, just a 
switch from one theoretical model to another, which can 
be justified, strictly speaking, only in the region where 
both models give adequate description of the system, i.e. 
the same result for the all possible observables. In this 
overlapping ST region such a transition can be realized 
at any infinitely narrow dividing hypersurface. On the 
other hand, for the case of a real physical phase transi- 
tion, like FO or/and hadronization, the infinitely narrow 
dividing hypersurface is an idealization of the layer of 
some finite thickness, as it was discussed above. We basi- 
cally assume that if we use the correct particle distribu- 
tions on the "post" side, we do not make a big mistake 
changing distributions to the new ones sharply instead of 
changing them gradually in some layer. 



